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PREFACE. 

This publication gives a general development of the 
theory of the Lambert conformal conic projection. It is 
intended to supplement the matter found in Special 
PubUcation No. 47 entitled, ''The Lambert Conformal 
Conic Projection with Two Standard Parallels." It is also 
supplementary in a way to Special PubUcation No. 49, 
which contains the Lambert projection tables for the region 
in France, and to Special Publication No. 52, which gives 
corresponding tables for the United States, since it gives 
as a whole the mathematical development of the theory 
upon which they depend. 

A short account of Lambert's life and work is given in 
the introductory paragraphs, followed by a few pages upon 
the subject of projections in general. 
2 
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GENERAL THEORY OF THE LAMBERT 
CONFORMAL CONIC PROJECTION. 



By Oscar S. Adams, 
Geodetic Computer, U. 8. Coast and Geodetic Survey, 



Since this publication is to treat of some of Lambert's 
work, it is altogether fitting that it should be prefaced with 
a short account of his life, especially with a statement of his 
significance in the domain of map projections. Johann 
Heinrich Lambert was born at Mulhausen in Alsace in 
1728. He was the son of a poor tailor and his education 
was entirely the product of his own exertions, due to a 
systematic course of reading. It was his regular custom to 
spend 17 hours pier day in study and writing. At the early 
age of 16 he discovered, in computations for the comet of 
1744, the so-called Lambert's theorem. During the latter 
part of his life he resided in Berlin, where he was much 
honored for his ability. It was in the application of 
mathematical analysis to the practical problems of life 
that he especially excelled. His untimely death occurred 
in 1777 in the forty-ninth year of his age. 

His contributions to mathematics were the series which 
bears his name, the conception of hyperbolic functions, his 
theorem in conies, and the demonstration of the incom- 
mensurability of X. Both Lagrange and Gauss used parts 
of his work L points of dep^ture for their investigations. 

Lambert's work in the field of map projections needs 
careful consideration. He was the first mathematician 
to make general investigations upon the subject of pro- 
jections. Those who preceded him in this work limited 
themselves to the development of a single method of pro- 
jection, principally the perspective, but Lambert con- 
sidered the problem of the representation of a sphere upon 
a plane from a higher standpoint and he stated certain 
general conditions that the representation was to fulfil, 

5 



6 U. S. COAST AND GEODETIC ST7BVEY. 

the most important of these being the preservation of 
angles or conformality and equal surface or equivalence. 
These two qualities^ of course, caoi not be obtained in the 
same projection. 

Although Lambert did not fully develop the theory of 
, these two methods of projection, yet he was the first to 
express clearly the ideas regarding them. The former, 
conformality, has become of the greatest importance to 
pure mathematics, but both of them are of exceeding im- 
portance to the cartographer. It is no more than just, 
therefore, to date the beginning of a new epoch in the 
science of map making from the appearance of Lambert's 
work. What he accomplished is of importance, not only 
for the generaUty of the ideas imderlying it, but also for 
his successful appHcation of them in methods of projection. 
The manner in which he attacks and solves any particular 
problem is very instructive. He has developed several 
methods of projection that are not only interesting but 
that are to-day in use among cartographers, the most im- 
portant of these being the conformal conic projection. 

The initial problem presented to us in map making is 
the representation upon a plane surface of the relative 
positions, sizes, and shapes of features that are found upon 
the curved surface of the earth. A perfect representation 
is impossible, since the surface of the earth being non- 
developable, can not be spread out in a plane. There are, 
however, many different ways of obtaining approximate 
representations, the theory and properties of which con- 
stitute the subject of map projections. 

The positions of the points upon the earth are usually 
defined by their latitude and longitude. Hence, if we can 
devise a suitable method of representing the meridians and 
parallels upon the sheet, the points can be plotted by their 
positions relative to these lines and the map can be con- 
structed. The term ''projection" is evidently used in a 
wider sense than that which is given to it in geometry. 
The majority of map projections are not projections in the 
geometrical sense — that is, perspective projections, orthog- 
onal projections, etc. — but merely a network of meridians 
and parallels that makes possible a one-to-one correspond- 
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en(5e between the places upon the earth and the points 
upon the map. 
Some of the thmgs to be desired in a map are as follows: 

1. Preservation of the shapes of the countries. 

2. Preservation of the relative sizes of the countries in 
their representation upon the map. 

3. The distance between places should be in constant 
ratio to their distances as indicated upon the map. 

4. A great circle upon the earth should be represented 
by a straight line upon the map. 

5. The latitude and longitude of any place should be 
readily found from its position upon the map. 

6. The ease with which a projection can be constructed is 
also to be considered from the practical standpoint. 

Only part of these things can be attained by any given 
method of projection. 

The scale of a map in any given direction is the ratio 
which a short distance measured upon the map bears to 
the corresponding distance upon the surface of the earth. 
The definition is Umited to short distances because the scale 
of a map generally varies from point to point. It would 
of course be desirable that the scale of the map should be 
correct in every direction at every point and constant for 
all parts of the map. This is impossible, however, since if 
it were true, the map would be a perfect representation of 
the spheroidal surface and could be fitted to it. In any 
given method of projection, therefore, some of the features 
to be desired must be sacrificed. 

The representation of the shapes of coimtries as nearly 
correct as possible is one of the most important fimctions 
of a map. A large coimtry can not, of course, be repre- 
sented without some distortion when considered as a whole, 
but small areas may be mapped by similar figures. A pro- 
jection that preserves the similarity of small areas is called 
orthomorphic. The representation of one surface upon 
another so as to preserve similarity of elements has been 
called by mathematicians conformal representation. An 
orthomorphic projection is therefore a conformal repre- 
sentation of the spheroidal surface of the earth upon a 
plane. Orthomorphic projections are in general not of 
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much use in map making imless the meridians are straight 
Unes. 

The Lambert conformal conic projection fulfills this re- 
quirement, and it has lately been brought into especial 
notice by the fact that it is the projection that is used for 
the battle maps in France. In this projection, as is the 
case with any conic projection, the intersecting cone gives 
a better map than is given by the tangent cone. In the in- 
tersecting cone we have a shortening of the scale between 
the standard parallels and a lengthening of the same at the 
top and bottom of the map. 

The following is a somewhat simplified development of 
the mathematical theory of the projection, including a 
determination of the necessary and sufficient conditions 
for a conformal mapping. Use is made of the tangent cone 
only in order to determine the two limiting cases of the 
projection. 

It is well known that a plane curve can be expressed 
in parametric form in such a way that x = (l>(t) and j/ = 
4^{t), <f> and ^ being some known fimctions of the variable t 
Thus the straight line passing through the point (a, h) with 
direction cosines a and p can be expressed in the form 

x^a+at,y = b + pt. 

A circle can be given in the form 

a; = a cos ty y = a sin t 

A space or skew curve can be similarly expressed as , 

x=*(0,j/=^(0,2=x(0. 

Such a curve may degenerate to a plane curve imder cer- 
tain conditions; the obvious one is that one of the func- 
tions is identically zero. 
A straight line in space becomes 

x = a + aty y=^h + pt, z = c+yt 

A circular helix is given by the equations 

x = a cos ty y=a sin ^, z=ht. 

When the coordinates are given in terms of functions 
of two variable parameters, the locus becomes a surface 
except under certain conditions that will be specified 
later. Thus, the tangent surface of a skew curve ex- 
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pressed by the parametric equations x==<l>it), j/=^(0, and 
2 = x(0 is given by the equations 

a;-0(O +s<l>'(t),y==4^{t) +s4^'{t), z=x(t) +sx'(0, 

the primes denoting differentiation with respect to t By 
employing s and t as independent variables, the surface, 
given is that swept out by a tangent to the curve as it 
moves along the curve. The surface of a sphere is given 
by the equations 

aj = asin^coss, y = a sin^sin^, BJidz = a cos<. 

To avoid some of the diflBicult questions connected with 
the theory of curves and surfaces, the functions that 
will be considered in this article will be limited to domains 
in which they are single valued, finite, contiauous, and 
differentiable. These limitations are permissible in the 
cases to be considered. Thus, in the equations of the 
sphere the variables are Umited to the domains 

In regard to the fimctions three cases may occur. In 
the first place all of the functions may be constant, thus 
giving a single point. In the second place, the fimctions 
may not be constant, but any two of them may be func- 
tions of the third; this would give a curve and not a surface 
since the three functions could then be expressed in terms 
of some single variable, Vy and the equations would become 

x=-X{v), 2/= Y(v)y and z=-Z{rj). 

As the third possibihty, at least two of the functions are 
independent of each other. In this case the locus repre- 
sented is a surface. Of course any two may be inde- 
pendent, but all that is necessary is that at least two bear 
this relation to each other. 

Let <^ and ^ be the ones that are thus independent. 
This requires that their Jacobian is not identically zero, or 

du' dv 
2)^ d^ 



du dv 



M=0 
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10 U. S. COAST AND GEODETIC SURVEY. 

It should be noted that a parametric representation of a 
surface is not unique. If such a representation has been 
found any number of other forms can be obtained by set- 
ting the u and v of the first form equal to two arbitrary- 
independent functions of two other variables 8 and L 
Thus the equations may be u = X(s, t) and v=ii{8, t) with 
the condition 

+ 



dx 


b\ 


ds' 


dt 


ds' 


dix 
dt 



If one of the variables is held constant while the other 
varies, a curve is traced out upon the surface. Thus there 
is a set of curves for t^ = constant, and another set for 
v = constant. These curves are called the parametric 
curves and u and v are called the curvilinear coordinates. 
By a change of coordinates as indicated above a set of 
curves can be determined such that the curves, u = con- 
stant, are everywhere perpendicular to the curves, v« 
constant. In the equations of the sphere, a =^ constant 
gives a meridian great circle and t = constant gives a circle 
of latitude or a parallel. In this case the curves are per- 
pendicular or orthogonal. In a more general way any 
functional relation between the variables u and v deter- 
mines a curve upon the surface. This can be expressed 
in the form of one coordinate as a function of the other, as 
v = g (u)f or in the form F (u, v) = 0. As a further con- 
sideration any ordinary differential equation of the first 
order between u and v determines a set of oo* curves upon 
the surface, since the integration of the same would give a 
functional relation between u and v containing one arbi- 
trary constant. 

The element of length upon a surface is given by the 

equation dS^ = dx^ + df+dz^. 

But ^^ ^^ 

du dv 
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hence 

, rt/dx dx .by by , bz b2\ , , 



+ 



m<^y<w]^-- 



To abbreviate the expression it is customary to use the 
symbols 

p^bx bx^.by by^,bz^ bz 
bu bv bu bv bu bv 

With these the equation for the element of length becomes 

dS^ = E du" + 2F dudv+G dv' 

If V is held constant, the element of length becomes 

d8^ = E du^, since dv — 

If the direction is taken positive in which u is increasing, 
we find 

dS^^^du 

Similarly for u = constant, dS^ = -^0 dv 

bx bi/ bz 

It via constant dx = ^du, dy = ^du, and dz = ^ du. 

The direction cosines of the tangent to the curve v = constant 

dx 
are given by the equations a = ^^, etc. 

Thesebecome, a = -7= r—f i5=-7= rr^y and 7 = -7= r — 

^ bu -y/E bu ^E bu 

Similarly the direction cosines of the tangent to the curve 

u= constant are given in the form 



a = 



1 bx 



bx ^, 1 by . , 1 bz 

^^' ^^m bi' ^""^ '^ ^in bi 
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If ^ is the anglQ between these two tangents or the angle 
between the curves, we have from the known relation, 
cos e = aa'+Pfi'+yy' 

1 / dx bx by^ dy dz bz \ 
""V^^^ dv"^du bv^bu bvj 

F 



It follows, then, that the parametric curves u and v 
upon a siuiace on which F is identically zero form an or- 
thogonal system unless the -^EG is identically zero at 
the same tune. Since when F is zero neither E nor G 
is zero for the usual real surfaces, it foUows that, when 

We shall now determine the conditions that are neces- 
sary and sufficient to permit a conformal representation 
of a surface upon a plane. One siuiace is said to be con- 
formally represented upon another when any infinitely 
small section of the first siuiace is represented by a similar 
figure upon the other. 

If a siKf ace is to be represented upon a plane according 
to some law and x, y, z are the rectangular coordinates 
of the siuiace and u^ v the rectangular coordinates of the 
plane, then for every pair of coordinates u, v ia the plane 
a set of coordinates x, y, z of the siKface must correspond; 
that is X, y, and z must be functions of u and v. 

We may then consider that the siKf ace with the para- 
meters u, V is to be represented on the plane with the 
rectangular coordinates u, v. The question then becomes 
imder what conditions such a representation is conformal. 
Let us suppose that (u, v), (u+du, v+dv) and (u + du, 
v + 8v) are three points upon the siuiace forming an in- 
finitesimal triangle. This will be represented upon the plane 
by the triangle (u, v), (u+du, v+dv) and (u + 8u, v + 8v). 
These two triangles must be similar but since their positions 
relative to each other do not come into consideration, 
no account is to be taken of signs either of sides or angles. 
Two things are then required, first that the sides of the 
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triangle extending from the point (u, v) of the siuiace 
must be in the same ratio as their representations In the 
plane extending from the point (u, v) of the plane, and 
secondly the angle included between the sides of the tri- 
angle upon the surface must be equal to the angle included 
by their images in the plane. That is, the square of the 
elements of length upon the surface must be proportional 
to the square of the elements in the plane and the cosine 
of the angle between the lines upon the surface must 
equal the cosine of the corresponding angle in the plane. 

But ds'=Edu^ + 2 Fdudv+Gdv'. 

ds^=E8u^+2 F8u8v+0 8v^. 

The quantities Ey F, and G would be the same in each case. 
They may be constants or fimctions of u and v but since 





Fig. 1.— Surface triangle and its representatida upon a plane. 

each curve is infinitesimal and starts from the same point 
(u, v)y they would be constants for this point. On the 
other hand; in the plane, we have 

ds^^du^-^-dv^ 
88^-=du^ + 8v^. 

The first condition is therefore that the proportion 

Edu^-h2 F du dv-\-G dv^ ^ du^-\-dv^ 
E 8u^+2 F du 8v + G 8v^^ du^ + 8v^ 

may be identically true. This evidently requires that 
E==G and F=0. Thus the parametric curves must be 
orthogonal. 

The cosine of the angle between the curves upon the 
surface must now be determined. The direction cosines 
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of the tangent at the pomt (u, v) of the curve dS upon the 
surface are determined by 

dx dy dz 
SS' dS' dS 

or dxj^dx-j 

ou ov 



etc. 



^Edu^+2 Fdudv + Gdv"' 

Also the du'ection cosmes of the tangent at the point (u, 
v) to the curve 8S upon the surface are in like manner 

bx^ , dx ^ by ^ , by ^ 

^E du^+2 F du ^+G 5v^' -yjE 5u^+2 F bu bv^-0 bv'' ® ^' 
Consequently, if a is the angle between the curves, we have 



cosa = 



•y/(Edu'+2Fdudv+Gdv^) (Ebu^+2F5u8v+Gbv^'' 
or 

E du du+ F (du 8v+8u dv) +G dv bv 



cos a= 



•yf(Edu^ + 2 Fdudv + Gdv^) (E8u'+2F8u8v+GSv*) 

In the plane the corresponding direction cosines are 

du dv 

-yjdu^+dv^' -y/du^ + dv* 

and 8u 8v 

If J. is the angle between these two lines, we have 

. du 8u+dv^ 

cos ^ = — . 

-y/idu' + dv') (du^+dv^) 

The second requirement is, therefore, that the equation 

E du gu4- F (du dv + du dv) +Gdv 8v 

'yl(Edu' + 2 F du dv+G dv") (E 8u^+2 F Su 8v + G dv^) 

_ du 8u+dv 8v 

" ^(du^ + dv") (M'+Sv^) 
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should be identically true. This requires as before that 
E=G, and F=0. Thus of the two conditions only one is 
necessary and one is also sufficient since each leads to the 
same results. It is thus seen that the linear element upon 
the surface must be such that it can be expressed in the 

If we have a surface with the linear element dS^ = 
Vdv? + Vdx^ in which F is a function of v alone and 
17 is a function of u alone, it can be expressed in 
the required form. Let us put it in the form d/S* = 

TJV I -jj- -^-y)' By changing the parameters, we may let 

B = I —F^ and 17 = | -7=* The linear element then becomes 

dS^^UVW^-drf). 

When the linear element is in this form, the surface is 
said to be expressed in terms of an isothermal orthogonal 
system of parameters and the net of u, v curves is said to 
form an isothermal orthogonal net. The surface is divided 
by them into an 00 ^ group of infinitesimal squares. 

After the surface has been expressed in this manner 
in terms of isothermal orthogonal coordinates, the general 
conformal representation of the surface upon a plane can 
be determmed. This at the same time determines the 
general conformal mapping of the u, v plane upon another 
plane with point to point correspondence. This general 
representation is given by the equation aj+iy=/ (u-{-iv), 

in which i denotes as usual V — l and/ (u-\-iv) is any ana- 
lytic fimction of the complex variable (u+iv). The 
element of length dS^ =rn?(du^ -^dv*) is represented in the 
plane by d^==d7?-\-dy^. Let the curve C be ^•epresented 
by the curve C^ in the plane. If d represents the angle 
which the curve makes with the curve v = constant, we 
will have dSt='mdu and 

dSv du 



cosd=^^= 



dS -^dv? +dv2 ' 
dv 



sin ^=-7--- — 3— 
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If 6^ is the angle that G^ makes with the axis of X, we have 

__ dx 

^ /».• ;« /I du+idv Idu-hidv 

€»*=»COS O + l sin 0=^ /7r:o-.-7 o =-%/- 



or 6***= 



■^Jdv^+dv^ \du—idv 
du+idv 



du—i dv 
Similarly, 

e^-cme -^isine^ dx-^idy ^ I dx-^-idy 
e -cose^.+isme^, ^ds^-^dy' ^dx-^idy 

-ift dx-\-idy 
or e^^= j 'J 
ax— I ay 

Therefore, 

du — idv dx-\-idy 

dx-\-i dy * du — idv 

But since x+iy=/(u+w), so also will x— ii/=/(u— iy) 

By differentiation, we have 

dx+idy=^ (du'\-idv)f^ (u+iv) 
or 

and 

dx—idy _.,, • x 

the primes denoting differentiation with respect to the 
complex variable. By substituting these values, we have 

If r and r^ are another pair of corresponding cmves 
starting from the same point and their angles are denoted 
by and 0„ we should find that 6^^*~^»> equaled the 
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same expression. Since this expression is constant for a 
given point, it follows that 

or 

6 — $1 = <l> — <l>i 
so that 

That is, the angle between the curves upon the surface is 
preserved in the representation on the plaue. If we had 
made use of the functional relation x+iy=F (u—iv) we 

should find the relation e2^<*+^) = pi/ L 4. ^^^ 
Therefore in the same way as above 

or 

That- is, the angles are equal but directed in the inverse 
sense, a difference that is not taken into accoimt in con- 
formal mapping. The most general conformal mapping 
is given then by the relations 

x+iy=f(u+iv) 
and 

x+iy^F(u—iv) 

The two other forms are merely conjugates of these, i. e., 

and 

x—iy^^'Fiu+iv). 

It is easily shown that these functional relations also 
preserve the proportionality of lengths. By differentiating 
the relation x+iy=f (u+iv) we find that dx+idy=^ 
(du+idv) f'{u+iv). By differentiating the conjugate of 
this, WQ ohtBXudx—idy =^ {du—idv) f {u—iv). 
When these results are multiplied, we obtain 

{dxY- (idyy^[(duy-(idvy]f(u+iv)f(u-iv) 
or 

da^'{-dy^==(du^+dv^)f'{u+iv)f(u—iv) 

But the element of length in the plane has the form 

d^=-dx^-hdy^ 

67341"— 18 8 
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Also on the surface dS^ = m* (du^ + dv^) 
By combining these relations, we find 

dS' 



or 



ds^^:^f{y'-^iv)f{u-iv) 



ds _ -^f (u + iv)f' (u -^^^0 
'3B~~ m 



The two functions /'(i^ 4- iv) and/'(u— iv) are conjugate 
complex functions and hence their product is real since 




Fio. 2.— Surface of revolution. 

both the sum and product of conjugate complex fimctions 
are real. The fact that the angle between the curves is 
held would in itself be sufficient to establish the comfor- 
mality of the representation, but we shall have need of 
the above expression which is called the magnification of 
the representation. 

If we have the equation of a plane curve given in the 
plane in the form a;=u, 3/ = 0, z=f(u), a surface of revo- 
lution with this curve as a meridian is given by the 
equations x=u cos v,y==u sinv, z =f(u) . For this surface 
E=l+[f(u)Y, F=0, and G'^u^, the prime denoting dif- 
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f erentiation with respect to u. The square of the element 
of length is therefore dS^= {l+\f'(u)f}du'+u'diP 



or 



We can let 



ds'^u\^±^g^du^+dv] 



w 



-J 



VTTLfM? 



u 



du 



Then dS^ = u^(dv)^+dv^). The surface is thus expressed in 
terms of an isothermal orthogonal set of coordinates. 

Another test for a conformal mapping is that the func- 
tions must satisfy the Cauchy-Riemann partial differential 
equations. We have 

x+iy=f(u+iv) 

^^f(u+iv) or §^=/(^+w) 






dx 



Hence dxdy 

With the form 
we obtain 



^ , dx by 

x+iy^Fiu—iv) 

i^-F'iu^iv). 



dv 



iF'(u--iv). 



Therefore 



dx 
du 



^and— = ^< 
dv dv du 
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These are the equations that detennine that a given 
complex function is a fimction of a complex variable in 
the accepted sense. The first set presersres the angle 
between two given arcs m magnitude and sign. The second 
set results in reversing the angle, or the original angle and 
the image have opposite signs. Both x and y are functions 
that satisfy Laplace's equation for two dimensions. Thus 
with the first set 

dt;2 bubv' 

But yy ^ d^y 

bvdu bubv 

Therefore by addition 

Also b^y , ^y_r. 

The same results may be found from the second set. 
From these results we may conclude that the most general 
conformal representation of one surface upon a plane is 
given by setting the complex variable in the plane equal 
to any analytic function of the complex variable formed 
from the isothermal orthogonal coordinates of the surface, 
or to any analytic fimction of the conjugate of this com- 
plex variable. In the first case direct equaUty of angles 
is foimd; in the second, the angles are equal but turned 
in opposite directions. 

If the surface to be represented is an ellipsoid of revolu- 
tion, the parametric equations may be chosen in the 
foUowmgform: 

a: = a cos Jf sin u, t/ = asin Jf sin u, z = h cos u. 

a is the semi-major axis, h is the semi-minor axis, Jf is 
the longitude and u the complement of the eccentric angle 
of the generating ellipse or the complement of the reduced 
latitude. 






» 
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The element of length upon the spheroid is given by the 
equation dS^ = dz^ -\-dy^ +d^. 
But 

dx^a cos M cos u du—a sin Jf sin u dM; 

dy^'a sia M cos u du +a cos Jf sin u dM; 

dz^ —i siaudu; 

hence the element of length becomes dS* = a* sin' u dW + 

{a? cos* u 4- &* sin' u) dv?. Thus £*== a' sin' t^, F= o, and 6 = 

a' — &' 
a' cos ' u + 6* sin' u. If « is put equal to €' the equation 

becomes 

dS^ = a' sin' u [dJ/ ' + (cot' t^ + 1 - €') dv?]. 

The equation of the generating elhpse is ^ven by 

a5 = a sin u; 
z^h cosu. 




Fio. 8.— Generating ellipse with angles X, p, and «. 

The tangent of the colatitude j>, is the cotangent of the 

angle which the normal makes with the X axis. The 

tangent of the angle which the tangent makes with the X 

. . dz 
axis 13^* 

The tangent of the angle which the normal makes with 
the X axis, or the tangent of the latitude L is equal to 
dx 

But 

— T— «tan L^-T cot u 



or 
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or 

tan L = cot 2> = t cot u 
hence 

- tanu==tan p 

a ^ 

or -yjl — t^ tan u = tan j>. 

, (1 — €2)cos'2> . , 1 

cos' u =~ — r — 2 > since cos* t^ =-r"r7 — 5" 
1 — €* cos' 2>' 1 + tan'w 

. , sin'i> 

sm' i^ = i — a a 
1— e'cos' 3> 

^ cos'u cos' J) 

1— «*cos*3) 

(cot* tt + 1 - «*)d«' = 7T — ^}~^T^t ■ » 
^ (1 — «* cos' p)* am* J) 

°^ 1 - «» COS* pL (1 - «" COS' p)» sin* pj 

Let 

(1 — €' cQs' 2>)sm 2> 
then 

fl— r__^£__l f g sin P <^y . ^ f — csin p dp 
""Jsinp 2J1 — €Cosp 2J l + €COS2>* 

By integration 

d=logtan|-|log(l-€Cosp)+Jlog (1 + 6 cos p)+logG^, 

If the Umits of integration are so chosen that G (the con- 
stant of integration) becomes equal to unity; we have 

r P /1+ccos j>\i"l 

or 

i 

T 



hence 



e^ 



tan 2. / 1+^cosy V 

^^2 \l^€COSpJ 
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The element of length now becomes 

1 — c^cos^pV / 

The parameters have been reduced to an isothermal 
orthogonal system. 

We can now determine any number of conf ormal repre- 
sentations of the spheroid upon the plane. All that is 
necessary is to make use of the relation 

x+iy=tf (M±ie) 

f being any analytic fxmction and either combination of 

signs being used. 

With the relation x+iy =/ ( Jf-i^),let/(v) = Ke^''. 

We now have 

j^ «ir+nog r tan I-. fJ±i52?P^T^ 



^Ee'-^l^i'i^^^r] 



xetiM 



= Ztan'£. (" i+^^^P yfcos IM+i sin Im). 
2 \1 — 6C0sp/ \ ^ 

By equating the real parts and the imaginary parts, we 
obtain 

2 \1 — € COS pj 

^ 2 VI — c cos pj 

In this projection the parallels become concentric cir- 
cles. The equation for the radius is 

r=VSM^=^tan'|- (\±1^^^. 
^ ^ 2 \1 — € cos p/ 

The meridians are represented by radii of these concentric 
circles. This method of projection is the one known as 
Lambert's conformal conic projection, first developed by 
John Heinrich Lambert, in his ' 'Beitr&ge zum Gebrauche 
der Mathematik," Berlin, 1772. It was later fully dis- 
cussed by Gauss. 
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If an angle z is assumed sucli that 



tan 1= tang, ('l±i^2?^^ 



£\2- 



€ coa p^ 

this angle will be very nearly equal to the complement of 
the geocentric latitude. 




FiO. 4.— Oeneratiiig ellipse with angles L and L\ 

From the equation of the generating ellipse (p. 21), the 
tangent of the geocentric latitude is given by 



But 



or 



hence 



-=- cot u=tani' 
X a 



tan i=T cotu 

cotu=- tanX. 
d 

tan i' = - cot u= -5 tan L, 



a and b being the semimajor axis and semiminor axis, 
respectively. 

Then to a sufficient degree of approximation z=^—L\ 

The value of z can be computed rigidly very conven- 
iently by assimiing an angle q such that cos 2=6 cos p. 
Then, since 



cot 



2 Vl- 



cos 



£ 



coaq 



tan^otan ? cot* ^ 



« <. .> 
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or 

log tan o ^^^? ^*^ 2 "^ * ^^S ^^* 2 * 

However, the approximate formula determines z to within 
a few tenths of a second. By using this auxiliary angle 
the equations become 

x^Z'tan' 2 ^^® ^^ 

y^Z'tan'o ^^ ^^ 

r=-K'tan' 2* 

With these values x is reckoned downward from the center 
of the concentric circles and y to the right of the central 
meridian if M is reckoned positive in that direction. 

K and Z are as yet left arbitrary constants. Z may be so 
determined that the ratio of the lengths of any two arcs of 
parallels on the map may be equal to the ratio of the 
lengths of the arcs that they represent. If iV is the radius 
of curvature perpendicular to the meridian, or the length 
of the normal prolonged to the minor axis, a radian of the 
parallel L^ has the length N^ cos L^\ in the same way the 
length of a radian of parallel L^ is N^ cos Z,. Conse- 
quently, the ratio of the lengths of these two arcs is repre- 
sented by 

J^i cos Lt 

iVa COSia' 

Since the A factor^ in the tables for geodetic positions is 

equal to 

1 



the ratio becomes 



iVsinl" 
A^ cos lix 



A^ cos ij 
The arc upon the map that represents the radian of parallel 
Li has the length lr^=lK tan' ^. The radian of parallel Z^ 

1 See United States Coast and Geodetic Sarvey Special Publication No. 8, entitled 
*'Formulae and Tables for the Computation of Geodetic Positions/' 
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is likewise represented hylr2^lK tan'^ -o** The ratio of 
lengths will be preserved if we have 

tan|/ ^^^^^' 

7 _^ log COS Lj —log cos La— log Ai H-log A^ 

z z 

log tan 2^— log tan -^ 

Kmsj now be determined so as to hold not merely the 
ratio of the arcs of parallels L^ and L, but also to hold the 
exact length of these parallels. This is an excellent method 
of determination for mapping such an area as that of the 
United States. In this way we should have 

ZZtan^ |i=»iV, cos L,^:^-^.f 
hence r t 

-rr cos L/* COS 1j» 

f^sB * as = 

A, sin 1'' Z tan' ^ ^3 sin 1'' Z tan' ^ 

The twofold determination serves as a check on the com- 
putation. 

With this determination of Z and K, we shall compute the 
expression for the magnification at any point. We 
employed as the form of the function/ 



dS 



m 



From the equation of the linear element on the ellipsoid 
on page 23 



m' 



VI — «' cos' p 



^_ Vl-e»C08»y 

asmp ' 
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But 

2 \1— € cosp/ 
therefore 

Kl tm^i 



fe = 



\1 — € COS p/ ^ ^ 



asinp 
or 

Pfi iVn cos Ln cos Z« 

r^ being the radius of the circle on the map that represents 
the parallel of Ln and p« being the radius of the parallel. 
The last form is obvious from the conditions. 

If the parallels to be held are chosen about ^ of the 
distance from the bottom and the top of the area to be 
mapped, the proper balance will be preserved. The 
upper and lower part of the map will then be about as 
much too large in scale as the central part is too small. 
The scale along L^ and L3 will be exactly correct. With 
this value of K one can tell how much any parallel is in 
error of scale by computing a radian of the parallel and 
the length of the arc which represents it on the map. This 
is just a statement of the equation 

Pn 

With this projection a map could be made of an area 
such as that of the United States so that it would not be 
in error of scale in any part of it by more than 1 1/5 per 
cent. A polyconic projection of the same area is in error 
of scale by as much as 6^ per cent in some parts. A 
Lambert projectioir for the United States to be evenly 
balanced should hold parallels 29° and 45°. The scale 
would then be just about 1 per cent short along the 37° 
parallel and 1 1/5 per cent long along the 49° parallel. 

In Special FubUcation No. 52 of the United States 
Coast and Geodetic Survey, are given tables of coordinates 
for a Lambert projection of the United States. The 
standard parallels w^re chosen as 33° and 45° to lessen 
the scale distortion of the middle section of the country. 
This scheme gives about 23^ per cent scale distortion 
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along the parallel of 25^ as the greatest on the whole map. 
Only a small part of Florida and of Texas are affected by 
this increase of scale error. 

The coordinates for mapping the parallels are most con- 
veniently computed using as origin the point where the 
parallel crosses the central meridian^ the central meridian 
Q bemg the F axis and a perpendicular to it 

the X axis. All the formulas required for 
computation are as follows: 

tan U = —^ tanZ 

, log COS X^~log cos Za~log ^^4-log A^ 

log tan -^ -log tan "2' 
cosXj cosi. 




Z= 



Ai sm 1^ I tan« -J A^ sin 1^ I 



tan' 



r=£'tan«| 
x=r sin IM 



y^r (1 — cosZlf) = 2rsin'-2 =x tan-^-' 



Fig. 7.— Beprasenta. 
tion of a small geo. 

If it is desired to make the computation 
considering the earth as a sphere, it is only necessary to 

let € = 0. z then equals the polar distance or 2=p = o — i; 

A =A= — ; — Ty- When these values are inserted in the 
* ^ a sm I'' 

above formulas, the correct forms are given for the sphere. 

The difference of the radii gives the spacing upon the 
central meridian. If the parallel at the top and the one 
at the bottom are constructed by determining the coordi- 
nates of their intersections mth the meridians, the me- 
ridians can then be drawn. They can then be subdivided 
as was done in the case of the central meridian. In this 
way the coordinates of the other parallels can be deter- 
mined without computation. 
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If the two parallels to be held approach each other indefi- 
nitely, we shall have to determine the limiting value of I. 
This can be done in a number of different ways, and 
in fact the computation was made in six different ways 
with the same result in each case. The value that was 
obtained is the same as that obtained by Gauss and cited 
by Forsyth in his "Theory of Functions of a Complex 
Variable." Especial pains were taken with this, because 
Germain in his "Traits des Projections" obtains an 
erroneous result. The value is here determined by the 
same method that Germain said that he used, not because 
it is the simplest but to illustrate the fact that the correct 
result can be so obtained. 

The problem is to so determine I that CD may have the 
same ratio to BD that AB has, or to determine the place 
where two successive parallels are held at the same ratio. 
This is the same as saying that the magnification shall be 
a minimum at this point and can be true only at this point. 

length of BD on the map is equal to dr 

length of BD on the earth = — "" ^ 3 

(1— c'cos* J))? 

length of CD on the map = Z(r+dr)dJf 

length of CD on the earth 

LVl~€^cos'p VVl — c^ cos' pjj 

[ g sin p a cos p dp ae^ sin' p cos pdp" ! 
. + . ~- 3 idM 

VI — €'COS'p Vl — ^'COS'j) (1 — e' COS'p)2" J 



=[ 



asmp a(l - €') cos pdp'l , , w 
. 3 + 1"" pM%. 

yl — €' cos' p (1 — €' cos' p) 2 J 



To meet the required conditions, we must have the propor- 
tion >^ ^"^ 

[ asinp a(l-€') cos pdp 'l ^ 
,.- . ^., yi-c'cos'p d-c' cos' p)lj 

dr '^ a(l - €^)dp 

(1 — c' cos' p)-2 
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Therefore, since 

"Xr (1 — e* cos' p) sin p' 
?™co3 p. 
The magnification is least at 

p = cos~" I or ati=Bin-' I. 
If this least value of the magnification is denoted by t, 
we have 

If Ha taken equal to cos p, and K is determined so as to 
hold the length of the parallel of colatitude p, we shall 
have the case of a tangent cone. 



Fio. B.— SlarcognphEO prof eotlon ol the nartheni tiBmispbvt. 
Qermain obtained the erroneous value 

J _ cos p — f* cos p(l — 3 sin' p) 

If I becomes equal to unity, we have po = o and the 
tangent cone becomes a tangent plane at the north pole. 
The equations now become 
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r=Ztaii|- 

x=r COS M 
y=rsm M. 

This gives a projection for the spheroid analogoxis to the 
st«reographic projection for the sphere. With the sphere 
the value of € is zero and z becomes the polar distance. 
We then have a perspective projection from the south 
pole upon the plane tangent at the north pole. 

If I becomes equal to zero the tangent cone becomes a 

cylinder tangent at the equator, and p© equals ~. To 

determine the values it becomes necessary to evaluate 
some limits. K becomes infinite but so that Kl is finite 
and ^—x finite: 



Z= 



Kl 



No cos Lq 

IteLD} ^ 
No cos Lq 

tan^f 



When 






then 

lim ZZ«iVo— a 

In the general formulas 



X- K tan' 2-. / l+«co8p \^ ^^ 
2 Vl — ccosj)/ 

y-Ztan«^. (\±±^^)hmlM. 
To evaluate the limit, let us write x in this form 

X- KCOB IMe^ 'o« [tan | • (1^^^^)^] 

or on developing the exponential 
.-2rcosZi/{l+Zlog[tan|. (l±i^)i]+0(P)} 
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0(P) denoting terms with P as the lowest power of I. But 

Xo^K 
hence 

xo-x- JC(1 -008 UO - ZZoos ZJf log[tan |.(1±1^)V] 

-KcmlMOiP). 

Taking the limit when po== 2' '^^ 

.log[«,^.(l^l^>]. 

Denoting Xo— x by x and substituting for p its value ^r—Lj 

we have 

x-olog,[tan(|+j). (l^l^)f] 

log« denoting the Naperian logarithm. 

y-intan'^. /l+€Cos^\^ sJnZ^^ 

Taking the limit when Po==^9 2==o 

M being of course expressed in arc. Interchanging x and 
y to give the coordinates as usually plotted, we have 

x^aM 

y=alog.[tan(|+J.). (1^1|^)t]. 

This is a projection of the spheroid analogous to the 
Mercator projection for the sphere. If c becomes zero we 
have the Mercator projection of the sphere. 

We thus find that the stereographic and the Mercator 
are special cases of the Lambert projection and are there- 
fore conformal. 

Certain points in a conformal projection may be singular 
points at which the conformality fails. This is the con- 
dition at the pole in the Lambert conformal conic projec- 
tion. The angles between the meridians are not pre- 
served, since the angle between two meridians is Zif instead 
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of If as Upon the eartli. At such points, if 'W=f(?) is the 
complex relation, -3- is equal to zero or to infinity, llie con- 
formality fails at the poles in a Mercator projection but is 
preserved at the central pole in a stereographic polar pro- 
jection. Since the angle at the center of the syston of 
circles in the Lambert projection is equal to IM, the 360° 
in longitude is mapped on a sector of a circle with central 
angle equal to 360° X I. Since 2 is in the usual cases lees 
than unity, the central angle will be less than 360°. If I 



Sia. fi.— Uareatoc projaotloa. 

is equal to ^, the angle of the sector would be 270°; for I 
equal to % the angle would be 240°. 

In the war zone in France, the maps are made in the fol- 
lowing manner. Instead of using the exact formulas that 
have been developed, approximations are employed. A 
cone tangent at 55 grades (49° 30') is first determined with 
the parallel of 55« as the central parallel of the map. 
Along the central meridian the parallels are spaced from 
the formula 
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P being the distance along the meridian on the earth and po 
being the mean radius of curvature at 55*^. This formula 
for Ar is a Taylor series development of Ar from the rigid 
formula for r, correct to the third power of fi} 

The radius for the parallel of 55* is taken as No cot 55*^, 
iVo being the radius of curvature perpendicular to the 
meridian. The values of Ar being added to or subtracted 
from this radius give the radii of the other parallels. AM 
sin 55*^ gives the arc along the parallel corresponding to 

l6flL5^ 




Fio. 10.— ^Graticule for the Lambert projection of the northern hemisphere with 2—^. 

the difference of longitude AM reckoned from the central 
meridian. 

After these values are computed the whole scale is re- 
duced by 1 part in 2037. This gives us approximately a 
cone intersecting at the parallels of 53«^ and 57*^. The 
whole map is then covered with a system of kilometer 
squares with origin at latitude 55«^ and longitude 6*^ east 

1 See United States Coast and Geodetic Survey Special Publication No. 47, p. 13. 



